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It is widely believed that at high Reynolds number (Re) all turbulent flows approach a state of 
"fully developed turbulence" defined by a unique, Re-independent statistics of the velocity fluctua- 
tions. Yet direct measurements of the velocity fluctuations have failed to yield clear-cut empirical 
evidence of the existence of fully developed turbulence. Here we relate the friction coefficient (/) of 
rough-conduit flows to the statistics of the velocity fluctuations. In light of experimental measure- 
ments of /, our results yield unequivocal evidence of the existence of fully developed turbulence. 



Turbulence is the unrest that spontaneously takes over 
a streamline flow when the flow is made sufficiently fast. 
A conspicuous manifestation of this unrest, and one that 
lends itself readily to theoretical analysis, is the advent of 
fluctuations in the velocity field of the flow. In a classic 
paper 0] , Kolmogorov made a few plausible assumptions 
to show that in flows of high mean velocity (or high Rc) 
the statistics of the velocity fluctuations should become 
asymptotically invariant to further increases in the mean 
velocity. The universal turbulent state defined by the 
limiting statistics of the velocity fluctuations was termed 
"fully developed turbulence," and its existence has be- 
come a widely held belief. Support for this belief has been 
sought in direct, hot-wire measurements of the velocity 
fluctuations; yet, as recent research has shown Q,y,BLb), 
the results remain inconclusive. In this Letter, we seek 
to prove the existence of fully developed turbulence by 
harnessing empirical evidence other than the direct mea- 
surement of the velocity fluctuations. We start with an 
outline of the intricacies of the problem. 

Kolmogorov studied the statistics of the velocity fluc- 
tuation at the lengthscale I, ui. He used dimensional 

(i.e., 

the mean value of (uA 



analysis to show that the structure functionj_M; 

^ 2 ) must take the form 



( Ul ) 2 =P[Re,l/L] (si) 2 / 3 , 



(1) 



where P is a dimensionless function of the dimension- 
less variables Re and l/L, Re 3> 1 is a Reynolds number 
of the flow, L is the largest lengthscale in the flow, e is 
the mean value of the rate of energy dissipation per unit 
mass, and I is confined to the inertial range, L ^> / ^> r), 
where r\ is the viscous (dissipation) lengthscale Q . Given 
that Re ^> 1 and l/L <C 1, it is natural to identify a 
plausible asymptotic scenario for Re — > oo and l/L — > 0. 
To that end, Kolmogorov assumed complete similarity 



lim 



p, 



with respect to Re and l/L, or limj> e ^ 
where p > is a constant prefactor pj. Under this as 
sumption, the leading term in the asymptotic expansion 



of (ui) 2 is (ui) 2 lt = p(el) /A (where the subscript "It" 
stands for leading term), which is independent of Rc, 
and may therefore be said to define a universal limiting 



state for Re ^> 1 — the state of fully developed turbu- 
lence . Shortly after the publication of Kolmogorov's 
paper in 1941, it was objected that the asymptotic sce- 
nario proposed by Kolmogorov, and customarily known 
as K41, could not account for intermittency — a phe- 
nomenon whereby the rate of energy dissipation per unit 
mass fluctuates around its mean value, e |9|]. To account 
for intermittency, Kolmogorov himself argued that the 
similarity with respect to l/L might be incomplete, and 
went on to assume the simplest type of incomplete simi- 
larity with respe ct to I / L Under this assumption, the 
leading term of [ui) 2 is (u/) 2 lt = p(sl) 2 / 3 U/L) a , where 
a is a constant intermittency exponent |ld| . This alter- 
native asymptotic scenario, known as K62, has led to a 
vast body of research on intermittency Q. Nevertheless, 
in K62 (w;) 2 it remains independent of Re, as was the case 
in K41, and therefore continues to define a universal lim- 
iting state for Re > 1. In this sense, K62 represents 
only a minor departure from K41. More recently, a ma- 
jor departure from K41 has been suggested on the basis 
of new experimental results 0, Q • These results indicate 
that even at very high Re the prefactor p is not constant, 
but subject to a discernible dependence on Re — a depen- 
dence that is marked enough as to cast doubts on the 
existence of fully developed turbulence. To analyze these 
results, Barenblatt and Goldenfeld 0] argued that there 
might be no similarity with respect to Re. Further, they 
argued that the form of the structure function at high 
Re should be invariant under a natural redefinition of Re 
and showed that for this principle of asymptotic co- 
variance to hold the structure function (and therefore p) 
must depend on Re only through In Re. Last, they wrote 
p(lnRe) =p +p 1 6 + o(6), where S = l/|lnRe| < 1 El 
Under these conditions, 



lut 



InRc 



(2) 



and the existence of a universal limiting state depends 
on the value of po > 0. If po > 0, then (ui) 2 lt — 
po (el) 2 / 3 (l/L) a , which is independent of Re, and a uni- 
versal limiting state exists, identical with that predicted 



by K41 (if a = 0) or K62 (i f a/0) 13]. On the other 
hand, if po = 0, then ( U/ ) 2 lt = Pl (el) 2 l 3 {l/ L) a / InRc, 
which depends on Re, and a universal limiting state does 
not exist. To decide between these alternative scenarios, 
Barenblatt and Goldenfeld computed best fits of J2J to 
hot-wire data from a large wind tunnel and the atmo- 
sphere 0, both for po > and po = 0, and concluded 
that the data were "not inconsistent with either of the 
two possibilities" 0] . In another attempt at deciding the 
matter, Sreenivasan [|| studied a large set of hot-wire 
data, and concluded that the prefactor p is "more or less 
universal, essentially independent of the flow as well as 
the Reynolds number" [1J]. Nevertheless, he noted that 
the scatter in the data was large, and that to evince the 
behavior of p at high Re one would have "to cover a wide 
range of Reynolds numbers in a single, well-controlled 
flow, and use instrumentation whose resolving power and 
quality remains equally good in the entire range;" unfor- 
tunately, "such experiments and efforts are not yet in the 
horizon at present" ||. To decide the matter, we intend 
to resort to experimental data on the friction coefficient 
of rough conduits, /. These data (e.g., [Lo^) appear to 
be well suited to our purpose: they contain very little 
scatter and show beyond doubt that, for any fixed wall 
roughness, the leading term of / is independent of Re at 
high Re. Further, the observed dependence of the lead- 
ing term of / on the roughness is described accurately 
by a well-known empirical formula. The only problem, 
to which we turn presently, is how to relate / to (u/) 2 . 

The friction coefficient of a conduit of circular cross- 
section may be defined as / = r/pV 2 , where r is the 
shear stress on the wall of the conduit, p the density of 
the liquid flowing through the conduit, and V the aver- 
age velocity of the flow. We seek to obtain an expression 
relating / to the structure function of 10) • Now J5J was 
originally derived under the assumptions of isotropy and 
homogeneity, but the turbulent flow in a conduit is both 
anisotropic and inhomogeneous. Nevertheless, recent re- 
search [JL6| nas established that applies as well to flows 
that are neither isotropic nor homogeneous. Further, if 
vi denotes the characteristic velocity of a turbulent eddy 

of size I, we may identify vi = \J (u;) 2 , where (it;) 2 is 
given by J2J with L = D (the diameter of the conduit, 
which sets the largest lengthscale in the flow) (0, chap. 
7). Thus, it follows from J5J that, regardless of the value 
of po, the smaller the eddy the lower its velocity. With 
these considerations in mind, we now seek to derive an 
expression for r, the shear stress on the rough wall of the 
conduit. 

Let us call S the wetted surface tangent to the peaks of 
the roughness elements of the wall, Fig. ^ (We assume 
roughness elements of uniform size r, as in Nikuradse's 
experiments ^H-) Then, for Re 3> 1, the shear stress is 
effected by momentum transfer across S. Above S, the 
velocity of the flow scales with V, and the fluid carries 
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FIG. 1: Immediate vicinity of a the roug h-conduit wall with 
roughness elements of uniform size r |25|. The dashed line is 
the trace of the wetted surface S. 

a high horizontal momentum per unit volume (~ pV). 
Below S, the velocity of the flow is negligible, and the 
fluid carries a negligible horizontal momentum per unit 
volume. Now consider an eddy that straddles the wet- 
ted surface S (one half of the eddy is above S, the other 
half below). This eddy transfers fluid of high horizontal 
momentum downwards across S, and fluid of negligible 
horizontal momentum upwards across S. The net rate of 
momentum transfer across S is set by the velocity nor- 
mal to S, which velocity is provided by the eddies that 
straddle S. Therefore, if v n denotes the velocity nor- 
mal to S provided by the dominant eddy that straddles 
S, then the shear stress effected by momentum transfer 
across S scales in the form r ~ p Vv n . Now the size of 
the largest eddy that straddles S scales with r, the size 
of the roughness elements. This eddy provides a veloc- 
ity v r = J (u r ) 2 normal to S, where (u r ) 2 is given by 
with I = r and L = D. Smaller eddies do provide a veloc- 
ity normal to S, but these velocities are overwhelmed by 
the velocity of the eddy of size r. Thus, v n ~ v r , and the 
dominant eddy that straddles S is the largest eddy that 
straddles S. We conclude that r ~ pv r V, and therefore 

/ ~ v r /v = ^T^/v. 

To complete our derivation, we relate e to V and D 
using the phcnomenological theory, which is based on 
two tenets pertaining to the steady production of tur- 
bulent (kinetic) energy: (1) The production occurs at 
the lengthscale of the largest eddies in the flow and (2) 
The rate of production is independent of the viscosity. 
From these tenets and the equality of production and 
dissipation, it follows that we can obtain a scaling ex- 
pression for e, the rate of dissipation of turbulent energy 
per unit mass of liquid, in terms of the velocity of the 
largest eddies (which ~ V) and of the size of the largest 
eddies (which ~ D) [l^. The largest eddies possess a 
kinetic energy per unit mass e ~ I^ 2 and a turnover time 
t ~ D/V. These eddies persist for a time t, whereupon 
they split into eddies of size ~ D/2, thereby transferring 
their energy to smaller lengthscales. For the steady state 
to be preserved, a new set of large eddies must be pro- 
duced at time intervals t, implying that e — e/t ~ V 3 /D 
[iH . By using s ~ V 3 /D, I = r, and L = D in and 

substituting the result in / ~ *J (u r ) 2 /V, we obtain 
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FIG. 2: A comparison of Nikuradse's data with Strickler's 
expression. The straight line is log(/) = — 0.54 + log(r/D)/3. 
The data points correspond to the highest Re (ps 10 6 ) tested 
by Nikuradse. 

an expression relating / to the parameters po, p\, and a 
of the asymptotic expansion of the structure function. 

Now consider a conduit of fixed roughness, r/D = 
const <C 1. If Po = the leading term of / depends 
on Re: / h ~ y^T (r/D) 1 / 3+Q / 2 (lnRe)- 1 / 2 . Thus, for 
Po = the friction coefficient vanishes asymptotically at 
high Re, a conclusion that is at odds with all experi- 
mental data on rough-conduit flows. On the other hand, 
if po > the leading term of / is independent of Re: 
fit ~ V^o (>/£>) 1/3+Q/2 - Thus, forpo > the friction co- 
efficient tends to a positive constant at high Re, a conclu- 
sion that is qualitatively consistent with all experimental 
data on rough-conduit flows. Further, in the case a — 0, 
fit ~ y/Po( r /D) 1 ^ 3 1 which we recognize as Strickler's em- 
pirical expression for the friction coefficient of a conduit 
of roue hness r/D at high Re [HIH. In the case a^O, 
fit ~ \fp~o~ (r/D) 1 ' 3+a ' 2 , which is a generalized form of 
Strickler's empirical expression that accounts for the ef- 
fect of intermittency. Given that the experimental data 
can be fitted very well even if a is set to zero [Fig. E] , we 
infer that 1/3 ^> |a|/2, or \a\ <§C 2/3, consistent with the 
available estimates of a |l2j . 

From the previous paragraph, we conclude that equa- 
tion J2J together with the experimental data embodied 
by Strickler's empirical expression allows us to establish 
that po > 0. Therefore, the leading term of the structure 
function (ui) 2 is independent of Re, and the universal 
limiting state known as fully developed turbulence exists 
at high Reynolds numbers. 

The logic of our reasoning so far has been the follow- 
ing. We have shown theoretically that in a rough con- 
duit limRo^oo / = unless Pq > 0. In other words, 
we have shown that po > is a necessary condition 
for the friction coefficient of a rough conduit to tend 
to a positive constant at high Re. Then, given the un- 
equivocal experimental evidence that in a rough conduit 
liniRe->oo / = const > 0, we have concluded that it must 
be that po > 0. Interestingly, po > is not in general a 



sufficient condition for the friction coefficient of a con- 
duit to tend to a positive constant at high Re: the conduit 
must be rough. To elucidate this statement, we recall our 
scaling expression for the shear stress, t ~ pv r V. This 
expression indicates that the momentum transfer is dom- 
inated by the eddies of size r — the same size as the rough- 
ness elements. Now in a rough-conduit flow of sufficiently 
high Re, r exceeds the viscous lengthscale, r > i). In fact, 
from r\ = ii 3 / 4 e~ 1 / 4 (where v is the kinematic viscosity) 
and e ~ V 3 / D we can write rj/D ~ [y/VDfl 4 = Re~ 3/4 
and conclude that for any given r the condition r > r\ 
holds at sufficiently high Re. If Re is increased further, 
7] lessens and newer, smaller eddies populate the flow 
and become jumbled with the older eddies. Yet the mo- 
mentum transfer continues to be dominated by eddies of 
size r, and / remains invariant. This argument explains 
the behavior of / for all rough conduits, no matter how 
small the the roughness elements. If r = 0, however, the 
condition r > n cannot be attained, even at extremely 
high Re. Thus in a smooth conduit |2l| the momentum 
transfer will always be dominated by the smaller eddies 
in the inertial range, whose size scales with rj\ since r\ 
lessens as Re increases, the momentum transfer will be 
dominated by ever smaller eddies as Re increases, and 
limRe^oo / = — even though po > 0, and the turbu- 
lence will be fully developed at high Re. To verify this 
conclusion mathematically, we study Q3J for r ~ rj. By 
substituting r]/D ~ Re -3 ^ 4 in place of r/D, we obtain 

In accord with our discussion above, Q indicates that in 
a smooth conduit the friction coefficient vanishes asymp- 
totically at high Re, whether po = or po > 0. Had we 
to decide between p = and po > on the basis of (@J 
and experimental data, the answer would not be clear- 
cut, because (In Re)" 1 / 2 varies but very slowly at high 
Re. Nevertheless, we have established previously that 
po > 0. It follows that /it - ^Re~ 1/4 - 3a/8 , which in 
the case a = coincides with Blasius's empirical expres- 
sion for the friction coefficient of a smooth conduit at 
high Re 13. In the case a ^ 0, f u - y^Re -174 " 3 " 78 is 
a generalized form of Blasius's empirical expression that 
accounts for the effect of intermittency. 

From the previous paragraph, we conclude that if po > 
(as we established before), our theoretical predictions 
are in accord with the experimental data embodied by 
Blasius's empirical expression. 

In a recent monograph 22] , Tsinober has drawn a dis- 
tinction between "possibly universal properties of small 
scale turbulence" and "quantitative universal properties 
of turbulence at high Re" 0] . The former were first envi- 
sioned in the early history of turbulence physics; they are 
embodied by the asymptotic scenarios of Kolmogorov, 
K41 and K62, that define the possibly universal limit- 
ing state of fully developed turbulence. The latter have 
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long been known to engineers; they are manifest in a few 
properties of turbulent flows that become independent 
of Re at high Re, including the drag of bluff bodies and 
the friction coefficient of rough conduits, among others 
23| . Here we have established a relation between the two 
24j. In particular, we have found that the existence of 
fully developed turbulence is a necessary condition for the 
friction coefficient of rough conduits to tend to positive 
constants at high Re, as seen in experiments. In addi- 
tion, we have found that the existence of fully developed 
turbulence is compatible with the experimental evidence 
that the friction coefficient of smooth conduits tends to 
zero at high Re. On the basis of our findings, fully de- 
veloped turbulence may be termed a "universal property 
of small scale turbulence," without qualification. 

We thank J. W. Phillips for kindly reading our 
manuscript and suggesting ways of improving it. 
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